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SUMMARY 

Demodulators  typically  include  low  pass  filters  to  recover  signals  from  wideband 
frequency  modulated  carriers.  Multispeed  analogue  tape  machines  normally  require  the 
carrier  frequency  to  be  set  to  a value  proportional  to  tape  speed.  Reproducing  hardware 
usually  includes  sets  of  filters,  with  bandwidth  proportional  to  tape  speed,  for  use  with 
the  demodulators.  By  dividing  the  frequency  of  the  reproduced  signal,  a filter  corresponding 
to  a lower  tape  speed  can  be  used  and  low  pass  filtering  of  recorded  data  signals  can  be 
readily  achieved.  However  some  unwanted  noise  components  can  be  generated.  These 
noise  components  are  examined  critically  and  a method  of  reducing  them  is  demonstrated. 


POSTAL  ADDRESS:  Chief  Superintendent,  Aeronautical  Reeeerch  Laboratories, 
Box  4331,  P.O.,  Melbourne,  Victoria,  3001,  Australia. 
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1.  INTRODUCTION 


« 


The  need  sometimes  arises  to  record  signals  over  large  bandwidth  using  frequency  modu- 
• lation  (FM)  techniques  and  subsequently  low  pass  filter  the  reproduced  signals.  Such  a require- 

ment has  recently  arisen  in  the  acquisition  of  vibration  data  from  an  aircraft  fitted  with  acceler- 
ation sensors.  There  are  two  main  areas  of  interest  relating  to  the  vibrations: 

(i)  To  obtain  information  on  the  total  vibration  environment  the  amplitudes  of  acceleration 
components  with  frequencies  up  to  about  2 kilohertz  (kHz)  are  required.  Time  correl- 
ation between  vibrations  sensed  at  different  points  is  not  a requirement  in  this  case. 

(ii)  Vibrations  with  frequency  components  up  to  about  30  Hz  are  of  interest  for  flutter 
analysis.  Amplitude  and  phase  correlation  between  vibrations  sensed  at  different  points, 
or  at  the  same  point  but  in  different  directions,  is  necessary. 

Normally  the  amplitude  plots  required  for  (i)  would  be  obtained  using  a spectrum  analyser 
and  the  analysis  required  in  (ii)  would  be  performed  with  the  aid  of  a digital  computer.  In  the 
latter  instance  the  analogue  vibration  data  must  be  regularly  sampled  and  converted  to  digital 
form.  To  prevent  the  introduction  of  aliasing  errors  in  this  case  low  pass  filtering  of  the  repro- 
duced data  is  necessary. 

As  tape  recording  tracks  are  usually  in  high  demand  in  applications  involving  the  acquisition 
of  vibration  data  on  aircraft,  it  is  expedient  to  use  the  same  sensor  to  provide  signals  to  meet 
both  requirements. 

By  dividing  the  frequency  of  the  reproduced  FM  signal  (before  demodulation)  the  data 
can  be  transferred  to  a lower  carrier  frequency  and  a standard  demodulator  filter  with  corres- 
pondingly lower  bandwidth  can  be  used.  Hardware1  for  such  frequency  division  has  been  manu- 
factured at  these  laboratories.  Because  the  demodulator  filters  are  usually  fairly  well  matched 
in  their  amplitude  and  phase  characteristics  their  use  for  data  signal  filtering,  with  require- 
ments such  as  detailed  in  (ii)  above,  represents  quite  an  advantage. 

This  report  indicates  how  the  standard  demodulator  filter  can  be  used  to  meet  the  filtering 
requirements  for  the  specific  case  cited  above.  However  the  techniques  described  can  be  readily 
extended  to  other  filtering  applications. 

When  the  demodulator  filters  are  used  in  this  non-standard  manner  for  data  filtering 
purposes  noise  components  of  significant  amplitude  can  be  generated  in  the  demodulator 
output.  A detailed  analysis  of  these  noise  components  together  with  details  on  techniques  which 
can  be  used  to  reduce  them  are  included  in  this  report. 


2.  DETAILED  EXAMINATION  OF  REQUIREMENT  FOR  LOW  PASS  FILTERING  AND 
TIME  SCALE  EXPANSION  OF  RECORDED  DATA 

Airborne  magnetic  tape  recording  equipment  and  ground  station  reproducing  equipment, 
available  at  these  laboratories  for  use  in  the  application  mentioned  above  (Sec.  I),  are  multispeed 
machines  fitted  with  Intermediate  Band2  FM  electronics.  Carrier  frequency  and  available 
bandwidth  are  given  as  a function  of  tape  speed  in  the  following  table. 


Tape  Speed  cm/s  (Approx.) 
ips 

IS2 

(60) 

76 

(30) 

38 

(15) 

19 

(7-5) 

9 5 
(3  75) 

4-75 
(1  875) 

Carrier  Frequency  kHz 

108 

54 

27 

13-5 

6-75 

3-375 

Data  Bandwidth  kHz 
{ ± 1 dB  (decibel)  passband) 

20 

10 

5 

2 5 

1-25 

0-625 

j 

i 


I 


I 

1 


i 


\ 


r1 

V 
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To  conserve  tape,  the  lowest  tape  speed  which  will  allow  adequate  data  bandwidth  to  be 
obtained,  is  normally  used.  To  accommodate  the  required  2 kHz  bandwidth  (Sec.  1)  a tape 
speed  of  19  cm/s  (7$  ips)  is  selected. 

For  the  flutter  analysis  under  consideration,  components  with  frequencies  up  to  about 
30  Hz  are  of  interest.  Low  pass  Altering  with  cut-off  frequency  above  30  Hz  (say  somewhere 
in  the  vicinity  of  60  Hz)  is  required. 

For  the  particular  computer  analysis  to  be  used  it  may  be  assumed  that  the  conversion 
rate  per  channel  is  limited  to  30  samples  per  second.  If  six  samples  per  cycle  of  the  maximum 
frequency  of  interest  (30  Hz  at  the  time  of  recording)  are  assumed  to  be  desirable  then  the  maxi- 
mum frequency  of  interest  must  not  appear  as  more  than  about  8 Hz  ( * 30/6  Hz)  at  the  input 
to  the  converter.  To  satisfy  this  latter  condition  it  is  essential  that  the  reproducing  time  scale 
be  expanded  by  at  least  a factor  of  4 relative  to  the  recording  time  scale. 


3.  USE  OF  FREQUENCY  DIVISION  TO  MAKE  POSSIBLE  LOW  PASS  FILTERING 
USING  STANDARD  MAGNETIC  TAPE  REPRODUCER  HARDWARE 

Intermediate  bandwidth  FM  systems  usually  employ  pulse  averaging*’4  demodulators  as 
these  are  very  suitable  for  handling  wideband  data.  Normally  a constant  duration,  constant 
height  pulse  is  generated  for  every  zero  crossing  of  the  reproduced  FM  signal.  Because  there 
are  two  zero  crossings  per  cycle  of  the  FM  signal,  the  carrier  frequency  is  effectively  doubled. 
The  train  of  pulses  thus  generated  is  taken  to  the  input  of  the  demodulator  low  pass  filter  which 
separates  out  a DC  component  (which  must  be  backed  ofT)  and  the  modulating  signal.  Various 
other  higher  frequency  sideband  components  are  present  in  the  input  to  the  filter.  In  most 
applications  separation  of  the  modulating  signal  from  the  unwanted  sideband  components  is 
the  chief  function  of  the  low  pass  filter. 

When  the  demodulator  filters  are  considered  for  use  in  the  filtering  of  recorded  data  their 
performance  outside  the  passband  takes  on  new  significance.  Because  the  frequency  of  modulating 
signal  components  may  extend  well  beyond  the  carrier  frequency  which  results  after  frequency 
division  the  generation  of  unwanted  sideband  components  with  frequencies  within  the  filter 
passband  can  take  place.  These  components  will  be  analysed  subsequently  (Sec.  3). 

The  eight  channel  frequency  divider1,  developed  to  facilitate  low  pass  filtering  of  recorded 
data,  provides  division  factors  in  binary  steps  from  I to  32  (i.e.  six  factors).  Division  by  a factor 
which  is  a power  of  2 was  adopted  simply  because  the  range  of  tape  speeds  available  for  the 
tape  machine  considered  follows  a binary  sequence.  In  general  division  by  factors  which  are 
not  a power  of  2 may  be  used  but  the  carrier  frequency,  which  results  after  division,  would  not 
then  correpond  to  standard  demodulator  settings.  In  this  latter  case  a DC  shift  with  restrictions 
also  on  dynamic  range  would  result  when  the  data  are  demodulated. 

It  is  convenient  to  define  the  “corner"  frequency  as  the  nominal  ± I dB  (decibel)  filter 
passband  limit  which  ranges  from  20  kHz  for  IS2  cm/s  FM  reproduction  to  623  Hz  for  4 -73  cm/s 
FM  reproduction  as  indicated  in  the  table  of  Section  2. 

To  see  whether  a filter  with  corner  frequency  in  the  vicinity  of  60  Hz  can  be  obtained  with 
the  aid  of  the  frequency  divider,  the  ratio  r/>  of  the  recording  band  limit  to  the  desired  filter 
corner  frequency  is  calculated. 

Thus  rD  = 2500/60  = 41-7  for  the  application  considered. 

However  the  particular  frequency  divider  available  will  allow  only  ratios  which  are  a 
power  of  2.  The  nearest  actual  ratio  rA  which  can  therefore  be  used  is 

rx  = 32  = 2* 

In  general  rA  = 2 V where  /V  is  the  frequency  division  switch  setting1. 

Using  a division  factor  of  32  the  low  pass  filter  comer  frequency  is  set  to  2300/32  Hz 
(i.e.  78  Hz). 

To  perform  the  filtering  in  one  pass  of  the  tape  without  the  introduction  of  an  intermediate 
“copy”  the  recorded  tape  must  be  reproduced  at  152  cm/s.  Division  of  the  FM  signal  frequency 
by  32  actually  transfers  the  data  to  a lower  carrier  frequency  (3373  Hz,  equal  to  108000/32  Hz). 
This  transfer  causes  high  frequency  components  to  be  irretrievably  lost  but  the  main  frequency 
selection,  with  any  associated  amplitude  and  phase  shifts,  occurs  when  the  signal  is  demodulated. 
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The  act  of  reproducing  the  data  at  152  cm/s  (60ips)  rather  than  19  cm/s  ips)  causes 
the  data  frequencies  to  appear  8 times  higher  than  they  were  at  the  time  of  recording.  Thus  the 
78  Hz  cut-off  frequency  appears  in  real  time  as  625  Hz (»  8 x 78  Hz)  and  the  30  Hz  maximum 
frequency  of  interest  as  240  Hz. 

If  the  analysis  system  could  handle  a data  bandwidth  of  625  Hz  it  would  be  sufficient  to 
reproduce  the  recorded  data,  perform  the  frequency  division  using  a factor  of  32,  demodulate 
using  a 4 -75  cm/s  (1{  ips)  filter,  and  digitize  the  data  without  producing  a filtered  tape  copy. 
In  that  case  8 seconds  of  recorded  data  would  be  digitized  in  1 second. 

For  the  flutter  analysis  the  digital  equipment  will  not  allow  a 625  Hz  bandwidth  to  be 
achieved.  As  indicated  in  Section  2 the  conversion  rate  is  limited  to  50  samples  per  second  and 
the  time  scale  at  the  time  of  the  digital  conversion  must  be  expanded  by  at  least  a factor  of  4 
relative  to  the  recording  time  scale.  If  a factor  of  4 is  adopted  the  figures  given  in  the  following 
table  apply. 


Frequency 
at  the  time 
of  recording 

Frequency  at  the 
time  of  the  digital 
conversion 

Samples  per  cycle  at  50  Hz 
sample  rate 

(at  time  of  digital  conversion) 

Signal  of  maximum 
frequency  of  interest 

30  Hz 

7 5 Hz 

6-7 

Signal  at  low  pass 
filter  corner  frequency 

78  Hz 

19  5 Hz 

2-5 

To  achieve  a time  scale  expansion  of  4 relative  to  that  at  the  time  of  recording  the  following 
operations  need  to  be  performed  : 

(i)  The  FM  signal  which  has  its  frequency  divided  by  32  for  filtering  purposes  needs  to 
be  demodulated  using  a standard  4 -75  cm/s  (1  j ips)  filter  for  intermediate  band  tape 
systems. 

(ii)  The  demodulated  output  from  (i)  needs  to  be  taken  to  an  FM  recording  module  aligned 
for  152  cm/s  (60  ips)  intermediate  band  (108  kHz)  recording  and  plugged  into  a tape 
copying  machine. 

(iii)  A “filtered”  tape  copy  is  to  be  produced  by  recording  on  the  copying  machine  at 
1 52  cm/s.  This  copying  process  sets  the  filtered  data  corner  frequency  to  A of  the 
available  bandwidth  at  this  recording  speed.  The  time  scale  “expansion"  relative  to 
that  at  the  time  of  recording  is  J. 

(iv)  The  filtered  tape  copy  is  to  be  reproduced  at  4 -75  cm/s  and  the  output  of  the  demo- 
dulator, aligned  for  that  reproducing  speed,  is  to  be  coupled  to  the  analogue  to  digital 
converter.  As  the  data  are  contained  in  the  lowest  & of  the  available  bandwidth 
negligible  changes  in  the  amplitude  or  phase  characteristics  should  arise  in  this  last 
demodulation  stage.  The  time  scale  expansion  relative  to  (iii)  is  32  and  that  relative 
to  the  time  of  recording  is  4. 

With  the  arrangement  detailed  above  it  takes  four  seconds  to  convert  every  one  second 
of  recorded  data  (referred  to  original  time  scale). 

Simultaneous  handling  of  multiple  channels  of  data  is  possible,  thus  allowing  time  cor- 
relation. 


4.  PERFORMANCE  OF  LOW  PASS  FILTER 

The  performance  of  a low  pass  filtering  system  meeting  the  requirements  detailed  in  the 
previous  section  has  been  evaluated  using  hardware  available  at  these  laboratories.  To  do  this 
a tape  recording  consisting  of  multiple  tracks  of  simultaneously  recorded  swept  sinewave  signals, 
and  a track  on  which  was  recorded  a voltage  proportional  to  frequency,  was  prepared.  Recording 
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using  FM  techniques  was  performed  with  a 7-track  model  FR1260  tape  system  manufactured 
by  Ampex  Corporation.  Full  amplitude  (producing  ±40%  carrier  deviation)  swept  sinewave 
signals  were  simultaneously  recorded  on  tracks  1,  3,  5 and  6 and  a voltage  proportional  to 
frequency  on  track  7. 

Using  the  hardware  arrangement  indicated  in  Figure  1 a filtered  tape  “copy”  was  produced 
at  1 52  cm/s  (60ips)  tape  speed  as  outlined  in  the  previous  section.  Information  reproduced  on 
tracks  3 and  7 was  copied1  without  demodulation  and  that  from  tracks  1,  S and  6 was  filtered 
prior  to  being  "copied”. 

Two  frequency  sweeps  were  recorded  on  the  original  tape  (and  subsequently  copied),  one 
extending  well  beyond  the  passband  limits  of  the  filter  and  one  just  encompassing  the  passband. 
When  the  filtered  data  were  reproduced  at  152  cm/s  (using  the  wide  frequency  range  sweep) 
together  with  a channel  of  data  which  had  not  been  filtered,  the  response  curves  of  Figure  2 
resulted.  Since,  for  the  particular  example  under  consideration,  the  passband  of  the  filter  had 
been  located  at  approximately  the  lower  & of  the  available  bandwidth  of  the  “copy”  channel 
the  relative  shape  of  the  filter  response  curves  (for  both  amplitude  and  phase)  would  be  virtually 
independent  of  the  speed  (and  hence  the  demodulator  filter)  used  to  reproduce  the  copy. 

A series  of  double  “humps”  (Fig.  2)  with  maximum  amplitudes  near  2 fc  (where  ft  is  the 
carrier  frequency  which  results  after  frequency  division  and  equals  3-375  kHz  for  152  cm/s 
reproduction)  is  to  be  noted  outside  the  passband  of  the  low  pass  filter.  These  “humps”  will  be 
considered  in  more  detail  in  Section  5. 

Using  the  reduced  frequency  range  sweep  (up  to  1 - 75  times  the  nominal  bandwidth  [response 
within  ± 1 decibel]  of  the  low  pass  filter)  the  low  pass  filter  response  of  Figure  3 was  auto- 
matically plotted.  Because  the  characteristics  of  the  various  demodulator  filters  are  not  perfectly 
matched  some  variation  in  response  curves  would  result  if  a different  filter  were  used. 

Using  the  signal  copied  (unfiltered)  on  track  3 (Fig.  1)  as  the  phase  reference,  the  phase 
of  the  filtered  signal  “copied"  on  track  5 was  plotted  as  a function  of  frequency  as  indicated  in 
the  lower  graph  of  Figure  3. 

The  extent  of  the  phase  mismatch  between  signals  passed  through  different  filters  is  of  great 
significance  especially  where  correlation  of  the  signals  is  important.  In  Figure  4 the  phase  mis- 
match between  signals  recorded  on  tracks  1 and  5 (same  head),  and  between  signals  recorded 
on  tracks  5 and  6 (different  head)  is  plotted.  For  the  specific  example  considered  here  the  filtered 
outputs  from  tracks  I,  5 and  6 match  to  within  2-5  degrees  for  0 to  60  Hz  bandwidth  where 
the  frequencies  are  referred  to  the  recording  time  scale  (19  cm/s  for  which  the  bandwidth  is 
2 5 kHz). 

5.  ANALYSIS  OF  UNWANTED  SIDEBAND  COMPONENTS 

As  indicated  in  Section  4 a series  of  double  “humps”  is  observed  in  the  amplitude  response 
(filter  output  versus  frequency  of  modulating  signal)  of  the  low  pass  filter.  Nodes  occur  at  2 /e, 
4 /,,  6 ft  etc.  where  fc  is  the  carrier  frequency  which  results  after  frequency  division  (3-375  kHz 
for  th.-  example  under  consideration).  In  each  case  antinodes  are  observed  at  a separation  /* 
from  the  node  frequencies.  For  each  antinode  fx  is  approximately  equal  to  0 2 /( (equal  to  about 
670  Hz  for  3 -375  kHz  carrier  frequency).  However  the  frequency  of  the  components  observed 
at  the  output  of  the  low  pass  filter  is  /*  in  each  case.  Since  fx  is  just  beyond  the  corner  frequency 
of  the  low  pass  filter,  components  at  this  frequency  effectively  fall  within  the  filter  passband  and 
filtering  of  these  unwanted  components  from  desired  components  is  impossible  after  frequency 
division  as  outlined  in  Section  4. 

In  the  low  pass  filtering  example  of  Section  4 a wideband  signal  was  used  to  modulate  a 
high  frequency  carrier  and  the  frequency  of  the  modulated  signal  was  later  divided  to  make 
the  carrier  frequency  compatible  with  a filter  having  lower  cut-off  frequency.  It  has  been  shown 
experimentally  that  identical  “humps”  occur  if  the  wideband  signal  L taken  directly  to  a modu- 
lator set  to  the  lower  carrier  frequency,  and  the  modulated  signal  is  later  demodulated. 

For  the  purpose  of  the  analysis  of  unwanted  sideband  components  it  is  therefore  adequate 
to  consider  a low  frequency  carrier  modulated  with  a wideband  signal  the  frequency  of  which 
extends  well  above  that  of  the  carrier.  When  such  a signal  is  taken  to  a pulse  averaging  demo- 
dulator (Sec.  3)  the  low  pass  filter  will  receive  at  its  input  a train  of  constant  width,  constant 
height  pulses  starting  at  each  zero  crossing  of  the  FM  signal. 
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In  order  to  examine  analytically  the  generation  of  sideband  noise  components  a Fourier 
analysis  may  be  performed  on  the  train  of  pulses  received  at  the  low  pass  filter  input.  Any 
sideband  components  (excluding  components  at  modulating  frequency)  with  frequencies  in  the 
low  pass  filter  passband  will  give  rise  to  unwanted  noise  components. 

Monson5  has  considered  the  frequency  spectrum  of  the  input  to  a low  pass  filter  used  in  a 
pulse  averaging  demodulator.  He  analyses  the  output  from  an  FM  multivibrator  system  and 
shows  that  there  is  an  equivalence  between  the  multivibrator  FM  and  pulse  position  modulation 
(PPM).  It  will  now  be  shown  that  such  equivalence  applies  equally  well  for  pulses  generated 
at  the  time  of  the  zero  crossings  of  a modulated  sinewave. 

To  a first  approximation  the  FM  output  from  the  tape  reproducing  head  may  be  considered 
to  be  a constant  amplitude  frequency  modulated  sinewave  signal.  In  practice  the  amplitude 
and  phase  will  be  frequency  dependent  but  the  net  effect  of  such  variations  on  the  time  of  the 
zero  crossings  will  be  small. 

The  instantaneous  radian  frequency  <o  of  the  modulated  carrier  signal  e(t)  may  be  ex- 
pressed as 

(i>  = (oc  -f  tod 

where  <oc  is  the  radian  frequency  of  the  unmodulated  carrier  and  w*  is  the  instantaneous  radian 
frequency  deviation. 

cud  may  be  expressed  as 

“»< t = A em (/) 

where  em(f)  is  the  modulating  signal  which  is  a function  of  time  / and  A is  a modulation  constant. 

The  instantaneous  phase  0 of  the  modulated  carrier  signal  e(i ) is  given  by 


Thus  we  have 

e(/)  = E sin  0 

= E sin  (wci  4-  f cod  dt  + ^o)  (1) 

where  E is  the  amplitude  of  the  modulated  carrier  signal  and  <f>o  is  a constant  of  integration 
dependent  on  the  initial  value  of  e(t). 

If  lr-i  and  l,  are  the  times  of  consecutive  zero  crossings  (two  per  cycle  for  the  demodulator 
under  consideration)  then 


where  Te  is  the  period  of  the  unmodulated  carrier  signal. 

For  PPM  the  value  of  each  instantaneous  sample  of  the  modulating  wave  is  used  to  vary 
the  position  in  time  of  a pulse  relative  to  its  unmodulated  time  of  occurrence.  If  sampling 
coincides  with  the  time  of  appearance  of  the  position  modulated  pulse  the  sampling  is  said  to 
be  “natural”.  For  PPM  with  natural  sampling  the  change  in  position  of  the  r**  pulse  is  given  by 

tr~tR  = [«  S(/r)l  To  (3) 


where  /«  is  the  time  of  occurrence  of  the  r‘h  pulse  when  there  is  no  modulation, 

S(tr)  is  the  value  of  the  modulating  signal  S(t)  at  time  tr, 
a is  a modulation  constant, 

and  T0  is  the  repetition  period  of  pulses  in  an  unmodulated  pulse  train. 

The  time  between  successive  pulses  is  given  by 

tr  - tr-l  = To  [1  + « {S(tr)  - S(/r_l)}]  (4) 

Comparing  equation  (4)  with  (2)  the  equivalence  between  the  low  pass  filter  input  and  a 
PPM  signal  is  readily  seen,  giving 


To  = 


Tc 

2 


For  the  present  analysis  a cosinusoidal  frequency  deviation  will  be  considered. 
Let  tot  = u>d  cos  wmt 

where  u>m  is  the  radian  modulation  frequency  and  too  is  the  peak  frequency  deviation. 


(5) 

(6) 


to  D 

In  this  case  a S(t)  = sin  <*W  (7) 

The  modulated  carrier  signal  e(t)  is  now  given  by 


e(t)  = E sin  ( u>ct  + 


<»d  . . 

— sin  <om/) 

“m 


(8) 


if  6 is  arbitrarily  set  equal  to  zero  at  t = 0 (i.e.  a zero  crossing  occurs  at  t = 0). 

The  spectrum  of  a PPM  signal  has  been  considered  by  a number  of  authors  including 
Krauss  and  Ordung,*  Black7  and  Carlson8.  As  the  sideband  components  of  interest  in  the 
present  investigation  are  generated  when  the  modulating  signal  frequency  is  either  in  the  vicinity 
of  'he  repetition  rate  of  the  unmodulated  pulse  train  or  in  excess  of  that  rate,  the  simplifying 
assumptions  which  are  commonly  used  are  not  valid  here.  The  analysis  by  Black  is  suitable 
and  will  be  considered  here.  However  a detailed  analysis  of  the  PPM  spectrum  is  not  given 
by  that  author,  but  an  expression  (apparently  derived  from  unpublished  work  by  Edson)  for 
computing  the  value  of  the  Fourier  coefficients  is  given.  Black  does  provide  a detailed  analysis 
of  the  pulse  duration  modulation  (PDM)  spectrum  for  natural  sampling  using  a double  Fourier 
series  approach  based  on  work  by  Bennett*  on  the  evaluation  of  modulation  products  which 
arise  when  two  signals  of  different  frequency  are  added  together  and  applied  to  a rectifier.  The 
close  analogy  between  PPM  and  PDM  is  well  known,  and  it  will  be  shown  that  the  spectrum 
of  the  former  can  be  readily  derived  from  that  of  the  latter.  Bell  and  Sergent10  also  consider 
the  spectrum  of  a PDM  signal. 

Consider  a PDM  signal  for  which  the  leading  edges  of  the  pulses  are  spaced  at  regular 
intervals  and  for  which  the  trailing  edges  are  varied  relative  to  the  leading  edges  according  to 
the  value  of  the  modulating  signal  at  the  time  of  the  trailing  edge  transition  (natural  sampling). 
Thus  the  trailing  edges  follow  exactly  the  same  relationship  (equation  3)  as  the  PPM  signal. 
The  following  expression  for  a PDM  signal  Fi(i)  with  nature?  «»mr!ing  and  trailing  edge  modu- 
lation is  given  by  Black. 


„ . . , M ® sin  mw't  » 

F\(t)  = k -f  - cos  uimt  + £ £ 

l m-l  mn  m- 1 


JtAmrrM) 


sin  {maict  — 2mrrk) 


<20 
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±®  Jn{mnM)  . ( , . 

V sin  [mulct  + nu>mt  — Imnk 

mtr  V 


mr\ 
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(9) 


where  2ir/o»e  (or  To)  is  the  repetition  period  of  the  pulses  (leading  edges),  2nk/a>e  is  the  duration 
of  unmodulated  pulses,  m and  n are  integers,  and  (M/2)  cos  u»mt  is  the  modulating  signal  as 
specified  in  the  following  expression  for  the  duration  Tr  of  the  rth  pulse. 
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= To 


, M \ 

k + — cos  wmt 


Jn{z)  is  the  Bessel  function11  of  the  first  kind  with  argument  z and  order  n.  In  integral11  form 


"Jo 


cos  (z  sin  <f>  — n<f> ) d<f> 


When  n is  an  integer  (which  is  the  only  case  we  need  consider  for  this  analysis) 

7_„(z)  = (— l)»y»(r)  = M-z) 

For  small  values  of  z as  given  by  Hildebrand12 

J»(z)  ~ 2^n\  Z*  <"> 

and  ■/-«(*)  ~ z*  (12) 

in  the  sense  that  the  ratio  of  the  two  quantities  connected  by  the  symbol  ~ approaches  unity 
as  z approaches  0. 

A leading  edge  of  one  of  the  PDM  pulses  is  assumed  to  occur  at  / = 0 and  the  pulse  height 
is  considered  to  be  unity. 

The  PDM  spectrum  of  equation  9 may  be  written  more  generally  as 

_ , . M ® sin  uict 

F\(t ) = k + — cos  wmt  + X 

2 «-=i  mn 

-Z  Z sin  | m (a,,/  - 2.*)  + n (13) 

where  the  first  term  is  the  DC  component,  the  second  term  is  the  modulation  component,  the 
third  term  represents  a sawtooth  wave  of  negative  slope  and  abrupt  transition  at  r = 0 (Krauss 
and  Ordung  consider  pulse  modulated  signals  of  various  forms  as  being  the  summation  of 
sawtooth  waves),  and  the  fourth  term  varies  with  the  modulation. 

Consider  now  a PPM  pulse  train  with  pulses  of  duration  2 7rk/(oc.  Assume  the  leading  edge 
is  modulated  according  to  natural  sampling  of  the  modulating  signal  (given  by  (M/2)  cos  u>mt  as 
before),  and  further  assume  a leading  edge  transition  occurs  at  t = 0.  The  time  function  Ft(t) 
of  the  PPM  signal  can  then  be  written  simply  as 

FtKO  = F\  ^ <uel,  wm  J , k^  — Fi  ^a>cr,  c< >mt,  o)  (14) 

where  Fi(wel,  wmt,  k)  = F i(t)  as  defined  by  equation  13. 

Substitution  into  equation  13  yields 


+ Z Z — 1 - — — 

I *-»  —00 


_ . . , , ( kn\  * ® 2Jn(mnM) 

Fi(t)  = k + M sin  — kn  sin  <■>«,  ( / ) + Z Z 

“e  \ “r/  m - 1 «--oo 

si"  {(”+"-)*’}'“  {—  ('  - £)  + ” [~  ('  - *')  - i]}  "5) 

If  the  modulating  signal  (M/2)  cos  <umt  is  changed  to  (M/2)  sin  o»mt  (equivalent  to  (M/2)  cos 
(cuMr  — */2))  the  modified  time  function,  Fa(r)  say,  is  obtained  by  replacing  wmt  in  equation  13 
with  <omt  — ir/2.  The  final  term  then  becomes  n {com(/  — kn/u>e)  — *}. 


F*t)  - (-!)•  I 


UnmM) 


The  DC  term  is  obtained  by  substituting  n = 0 and  letting  m approach  0 in  equation  16, 
and  using  relationships  11  and  12. 

The  spectrum  of  the  low  pass  filter  input  can  now  be  easily  written  down.  From  equations 
3,  7 and  10  the  shift  of  the  rth  pulse  is  given  by 

oidT 0 , 

lr  — tR= sin  a >mt 

u>dT0  ( w\ 

= COS  ' a >mt  + - } 
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Thus  the  required  spectrum  F«(r)  is  obtained  by  substituting  2<oc  for  a*e,  2/3/w  for  M and 
<umt  + w/2  for  <«mr  in  equation  I S. 

+ -£)}  + £ £ 

w 2<oe  ( \ 2ci*c/J  *-l  ■ Wf 

sin  {(” + * £) } cos  { *"“■  (■  -£)+”"■('  - £)}  <|7) 

If  the  zero  time  reference  (r  = 0)  were  shifted  to  the  middle  of  the  first  pulse  then  t — kn/a>c 
in  equation  17  would  be  replaced  with  l. 

The  ( — 1)"  factor  in  the  PPM  spectrum  of  equation  16  has  not  carried  through  to  the  FM 
demodulator  low  pass  filter  spectrum  of  equation  17  and  indeed  would  not  have  been  included 
in  equation  16  if  the  modulating  signal  had  been  defined  as  — (M/2)  sin  wmt.  Although  it  was 
of  no  consequence  in  his  analysis,  Monson  seems  to  have  overlooked  this  point  as  he  indicated 
a negative  sign  for  the  component  at  the  modulating  frequency. 

For  convenience  equation  17  may  be  replaced  by  the  following  equations 


F«(r)  = + ^oi  cos 


M-?)} 
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where  ft  = u>e/2ir,  fm  = u>m/2 n,  fD  = wdI2w  and  to  is  the  pulse  duration.  Note  that  now 
2 trie  = (2<uc)fo  or  k = 2 fc  to- 
ll was  indicated,  for  the  filter  example  under  consideration,  that  for  the  major  antinodes 
in  the  response  (filter  output  versus  modulating  signal  frequency)  the  observed  filter  output 
frequency  was  /*  (approximately  equal  to  0-2 /c  or  670  Hz).  The  various  values  of  fm  required 
to  yield  a significant  sideband  component  at  fx  frequency  and  the  relative  amplitude  of  such 
components  are  of  interest  and  will  be  tabulated  shortly. 

The  required  sideband  frequency  can  be  expressed  as 

fx  = I 2m ft  4-  nfm  I (20) 

and  values  of  modulating  frequency  which  produce  that  sideband  frequency  are  given  by 


f 


The  normalized  amplitude  of  the  relevant  sideband  components  as  given  in  equation  23 
is  tabulated  below. 


* A negative  sign  simply  indicates  the  initial  phase  of  the  cosinusoidal  spectral  component 
(equation  17). 


The  analysis  of  this  section,  the  results  of  which  are  summarized  in  the  preceding  table, 
predicts  that  the  most  significant  sideband  component  will  be  generated  when  the  modulation 
frequency  is  1-8 fc  (equal  to  6 08  kHz  for  3-375  kHz  carrier  frequency)  and  that  it  will  have  a 
relative  amplitude  of  about  10-8%.  The  calculated  amplitude  of  the  component  which  arises  at 
this  modulation  frequency  and  of  those  at  other  antinodal  frequencies  agrees  closely  with  the 


measured  response  as  detailed  previously  in  Figure  4.  Some  slight  reduction  in  the  amplitude 
of  the  fx  sideband  component  is  to  be  expected  as  it  passes  through  the  low  pass  filter  since 
the  response  has  started  to  “roll  off"  at  frequency  /*. 

Consider  now  /*  as  a variable  frequency  as  defined  by  equation  20  for  n = — 1 . When  fx 
is  in  the  vicinity  of  the  filter  corner  frequency  (e  g.  ft  = 0 2 /<•  as  considered  above)  the  filter 
output  will  appear  sinusoidal.  When  /*  is  much  lower  than  the  filter  comer  frequency,  the  filter 
output  will  appear  non-sinusoidal  (noted  experimentally)  because  of  the  presence  of  harmonics. 
As  an  example  the  summation  of  the  particular  spectral  components  listed  in  the  following 
table  results  in  the  generation  of  a non-sinusoidal  signal  of  fundamental  frequency  /*  at  the 
output  of  the  filter. 


Value  of  m 

Value  of  n 

Frequency  of 
Spectral  Component 

Comment 

1 

-1 

/«.  2fc  -/. 

Fundamental 

2 

-2 

2f„  4/e  - lfrn 

Second  harmonic 

3 

-3 

3/„  6/e  - 3 fm 

Third  harmonic 

As  ft  -*  0 (i.e.  fm  -*■  2 mfe)  the  amplitude  of  the  spectral  components  (as  predicted  by 
equation  22)  at  frequency  /,,  and  of  harmonics  of /*,  approach  zero  thus  giving  rise  to  the  nodes 
noted  in  Figure  2.  The  transmission  through  the  filter  will  increase  with  separation  of  the  modu- 
lating signal  frequency  fm  from  the  nodal  modulating  frequencies  2mfc  until  limited  by  attenu- 
ation outside  the  passband  of  the  filter.  Thus  the  filter  output  corresponding  to  the  antinodes 


10 


of  Figure  2 is  expected  to  exhibit  a frequency  equal  to  that  at  which  the  / A(J)  function  (where 
A(f)  is  the  amplitude  response  characteristic  of  the  filter)  peaks.  The  observed  0-2  fe  (670  Hz 
approximately  for  the  particular  example  considered)  agrees  with  this  prediction. 


& METHOD  OF  REDUCING  SIDEBAND  NOISE 

Excessive  data  bandwidth  at  the  time  of  the  original  recording  may  give  rise  to  an  un- 
acceptable level  of  unwanted  sideband  components  if: 

(i)  The  modulating  signal  frequency  at  the  time  of  the  original  recording  is  greater  than 
about  1 - 6 times  the  carrier  frequency  (equivalent  to  the  leading  edge  of  the  first  sideband 
“hump"  in  Figure  2). 

(ii)  The  modulating  signal  frequency  at  the  time  of  the  original  recording  is  greater  than 
about  1 - 6 times  the  carrier  frequency  arising  after  frequency  division  (which  precedes 
low  pass  filtering  of  the  type  described  in  Sec.  4). 

Consider,  as  an  example  of  (i),  an  application  in  which  a wideband  transducer  is  used  for 
vibration  measurement,  but  only  the  lower  frequency  components  are  of  interest  for  subsequent 
analysis.  Assume  that  a tape  speed  and  hence  also  carrier  frequency  are  chosen  such  that  the 
corresponding  demodulator  filter  should  restrict  the  data  passband  to  the  band  of  interest  only. 
As  a specific  example  assume  that  components  with  frequencies  in  excess  of  10  kHz  are  present 
in  the  transducer  output  at  the  time  of  recording.  Further  assume  that  a 625  Hz  cut-off  frequency 
(as  obtainable  with  3-375  kHz  carrier  for  4 -75  cm/s  recording)  adequately  encompasses  the 
frequency  band  of  interest.  A vibration  component  at  6-08  kHz  frequency  would  give  rise  to 
an  unwanted  670  Hz  component  at  about  10%  the  amplitude  of  the  6-08  kHz  component  at 
the  time  of  recording.  In  cases  where  the  level  of  the  unwanted  sideband  components  which  pass 
through  the  reproducing  filter  is  excessive,  frequency  band  limiting  prior  to  recording  would  be 
essential. 

A similar  problem  to  that  detailed  in  the  previous  paragraph  may  arise  if  recording  is 
performed  with  an  adequately  high  carrier  frequency  but  frequency  division  and  low  pass 
filtering  are  performed  when  the  data  are  reproduced  (an  example  of  (ii)).  However  in  this 
latter  case  a modified  frequency  division  and  filtering  procedure  can  be  employed  to  reduce 
the  level  of  the  unwanted  sideband  components  if  that  is  considered  necessary. 

The  first  major  sideband  peak  occurs  when  the  modulating  signal  frequency  is  about 
6-08  kHz,  if  a straight  division  by  a factor  of  32  is  employed  for  the  example  considered  in 
Section  4.  For  the  modulating  signal  frequency  below  about  5 -4  kHz  the  sideband  contribution 
is  virtually  zero.  It  follows  therefore  that  if  the  upper  modulating  signal  frequency  is  kept  below 
about  5-4  kHz  (equal  to  1-6  times  carrier  frequency  or  8-6  times  filter  comer  frequency  for 
4-75  cm/s  reproduction)  the  sideband  components  should  be  largely  eliminated.  The  latter 
condition  can  be  satisfied  if  the  data  are  pre-filtered  prior  to  reaching  the  final  filter.  In  Figure  5a 
the  hardware  arrangement  for  low  pass  filtering  as  discussed  earlier  is  drawn,  and  in  Figure  5 b 
a modified  arrangement  which  provides  pre-filtering  with  2-5  kHz  corner  frequency  is  drawn. 
An  additional  stage  of  modulation  and  demodulation  (13  5kHz  carrier  frequency)  using  standard 
hardware  is  required  following  a frequency  divider  with  division  factor  of  8.  The  final  frequency 
division  factor  is  4. 

A comparison  of  the  amplitude  response  curves  obtained  experimentally  when : 

(i)  Pre-filtering  is  not  included 

(ii)  Pre-filtering  is  included 
is  shown  in  Figures  6 and  7. 

Reduction  of  sideband  components  with  pre-filtering  is  readily  observed  in  Figure  6. 
About  2%  maximum  sideband  level  is  noted  for  a modulating  signal  frequency  near  27  kHz. 
Significant  sideband  components  will  be  generated  at  the  1 9 cm/s  (7J  ips)  pre-filter  output  in 
the  vicinity  of  27  kHz  modulating  frequency  (twice  carrier  fluency  for  19  cm/s  recording). 
The  amplitude  of  these  sideband  components,  generated  as  the  modulating  signal  frequency  is 
varied,  will  increase  from  zero  at  zero  sideband  frequency  to  peak  at  about  2-7  kHz  (0-2/e) 
sideband  frequency.  Components  with  frequencies  falling  within  the  passband  of  the  4-75  cm/s 
(I  | ips)  filter  will  give  rise  to  the  sideband  components  noted  in  Figure  6 near  27  kHz  modulating 
frequency. 


II 


Some  modification  to  the  low  pass  filter  amplitude  response  (Fig.  7)  arises  as  a result  of 
the  additional  stage  of  pre-filtering  but  for  the  particular  filters  used  the  difference  is  very  small. 
Some  modification  to  the  phase  response  will  also  result. 

To  reduce  the  level  of  unwanted  sideband  components  and  also  to  avoid  excessive  departure 
in  the  low  pass  filter  characteristic  from  the  normal  response  a final  frequency  division  factor 
of  4 (even  if  a division  factor  other  than  32  is  required)  represents  the  optimum  value. 

7.  CONCLUSION 

a)  Low  pass  filtering  of  data  recorded  on  magnetic  tape  using  frequency  modulation 
techniques  can  be  simply  performed  when  the  data  are  reproduced  using  a frequency 
divider  and  standard  demodulator  filter. 

b)  The  filtering  technique  described  can  give  rise  to  the  generation  of  unwanted  sideband 
components  having  frequencies  within  the  filter  passband. 

c)  Analytical  expressions  for  the  level  of  unwanted  sideband  components  can  be  derived 
from  a Fourier  analysis  of  the  filter  input  signal. 

d)  A pre-filtering  technique  can  be  used  to  greatly  reduce  the  level  of  the  unwanted  sideband 
components. 
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FIG.  2 COMPARISON  OF  FILTERED  AND  UN-FILTERED  DATA  PASSBANDS 


FIG.  3 AMPLITUDE  AND  PHASE  RESPONSES  OF  LOW  PASS  FILTER 


FIG.  4 EXTENT  OF  PHASE  MISMATCH  BETWEEN  VARIOUS  FILTERS 


FIG.  5 (b)  HARDWARE  WHICH  INCLUDES  PRE-FILTER  TO  REDUCE  SIDEBAND  NOISE 
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FIG.  7 EFFECT  OF  PRE-FILTERING  ON  LOW  PASS  FILTER  AMPLITUDE  RESPONSE 
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